Role of d orbitals in the Rashba-type spin splitting for noble-metal surfaces 
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We investigate the Rashba-type spin splitting in tlie Sliockley surface states on Au(lll) and 
Ag(lll) surfaces, based on first-principles calculations. By turning on and off spin-orbit interaction 
(SOI) partly, we show that although the surface states are mainly of p-orbital character with only 
small d-orbital one, d-channel SOI determines the splitting and the spin direction while p-channel 
SOI has minor and negative effects. The small d-orbital character of the surface states, present even 
without SOI, varies linearly with the crystal momentum k, resulting in the linear k dependence of 
the splitting, the Hallmark of the Rashba type. As a way to perturb the d-orbital character of the 
surface states, we discuss efi'ects of electron and hole doping to the Au(lll) surface. 
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I. INTRODUCTION 

Approximated to the second order in power of the in- 
verse of the speed of light, c, the relativistic Dirac equa- 
tion is reduced to the usual nonrelativistic Schrodinger 
equation for the two-component spinor wave function 
with the three relativistic correction terms: the mass- 
velocity term, the Darwin term, and the spin-orbit in- 
teraction (SOI) term.— Among them, SOI induces many 
intriguing magnetic and electronic phenomena in con- 
densed matter physics, such as the intrinsic spin Hall 
effeclj^i^ and the Rashba-Bychkov effect (simply, the 
Rashba effect) in two-dimensional (2D) semiconductor 
heterostructures,— and recently it is known as one of the 
main mechanisms of the generation of the robust metallic 
surface (or edge) states in topological insulatorsi^ These 
phenomena have been expected to open a new pathway 
to the spin-dependent electronics and the quantum com- 
puting, and thus SOI in condensed matters has been of 
great interest both theoretically and experimentallyi^ 

So far, exploration of SOI phenomena has been lim- 
ited to materials with elements of high atomic numbers 
such as noble metals of Pt and Au, and p-orbital (scmi- 
)metals like Sb, Tl, Pb, and Bi. Among them, Au has 
the atomic SOI strength, (1/2)^5^, of 0.35 eV for the 
5d level,— >^ attracting special interest due to the SOI- 
induced Rashba-type spin-splitting of the surface states 
in Au(lll) surface first observed by angle-resolved pho- 
toemission spectroscopy.— Such spin splitting was also 
observed in other clean (semi-)metallic surfaceai^^ii and 
surface alloysj ^^'^'^ however, the splitting was unresolv- 
ably small in another noble-metal Ag surfac o^"*^'^^ al- 
though Ag has the atomic SOI strength, (l/2)^4c;, of 
0.13 eV for the 4d level.i 

Previous theoretical studies pointed out that the dis- 
similarity in Au and Ag is attributed to their differ- 
ences in the intrinsic atomic SOI strengtbi^ and the de- 
gree of asymmetry of the wave function near the nucleus 
of surface atoms or, stated differently, the hybridiza- 
tion of different-parity orbitals in the surface-state wave 
functions j^ii^ Compared with Ag, Au has a distinguish- 
ing feature that the energy-level separation between Sid 



and 6s orbitals is rather small,—"— which is related to 
the relativistic mass-velocity and Darwin termsJ^ This 
proximity of 5d- and 6s-orbital energy levels in Au in- 
creases d-orbital character in the surface states,— which 
is found important for the splitting.— However, the 
role of d orbitals in the splitting is still unknown except 
that they make the surface states asymmetric. 

In our present work, we investigate the origin of the 
Rashba-type splitting in the Shockley surface states on 
Au(lll) and Ag(lll) surfaces, focusing on the role of d 
orbitals. By turning on and off parts of the SOI Hamilto- 
nian, we find that the splitting originates from interaction 
between d orbitals via d-channel SOI although the surface 
states are derived mainly from p orbitals with only small 
d-orbital character, and the spin direction is also deter- 
mined by the d-channel SOI. The linear dependence of 
the splitting on the wave vector, k, which is the Hallmark 
of the Rashba type, is found to be derived from fc-linear 
dependence of dxz (dj,z)-orbital character of the surface 
states which exists even without SOI. These results show 
that the small d-orbital character of the noble-metal sur- 
face states determines the Rashba-type splitting while 
the large p-orbital one has weak and negative effects. As 
a way to modify the splitting by varying d-orbital char- 
acter of the surface states, we consider electron and hole 
doping to the Au(lll) surface. 



II. CALCULATIONAL METHOD 

Our present calculations are based on ab initio density- 
functional methods)^!^ which employ a Troullier- 
Martins norm-conserving pseudopotential^^ and the 
Perdew-Burke-Ernzerhof-type generalized-gradient ap- 
proximation for the exchange-correlation potentiali^ 
Wave functions are expanded with the localized 
pseudoatomic orbital bases.— The SOI term is in- 
corporated within a fully relativistic j-dependent 
pseudopotential)^ and treated in the ^-dependent fully- 
separable nonlocal form using additional Kleinman- 
Bylander projectorsi^^^ In treating the SOI Hamilto- 
nian, we consider non-self-consistently only the effect of 
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atomic core potentials with the L ■ S form, where L and 
S are the orbital and spin angular momentum operators, 
respectively. Although we can also include contribution 
of self-consistent part of potentials to SOI, our test cal- 
culations confirm that this contribution makes very lit- 
tle difference in the electronic band structuresi^ which is 
smaller than 10~^ eV. Thus, we do not include the con- 
tribution of the self-consistent part of the potentials to 
SOI in our following calculations. 

For surface calculations, we model Au(lll) and 
Ag(lll) surfaces with supercell including a slab of 30 
atomic layers and a vacuum region of about 20 A thick- 
ness, and we take 36x36x1 Monkhorst-Pack special- fc- 
point meshes in the full Brillouin zone (BZ)^ to inte- 
grate the charge density. While Ag(lll) exhibits no sur- 
face reconstruction, the Au(lll) surface reportedly has 
a surface reconstruction of a 22x-\/3 unit cell, the "her- 
ringbone reconstruction" ^ Since this reconstruction has 
negligible effect on the electronic structures of the surface 
stateS)2£ we safely take the clean-cut bulk atomic struc- 
tures for Au(lll) and Ag(lll) surfaces<^ 



III. ELECTRONIC BAND STRUCTURES FOR 
Au(lll) AND Ag(lll) SURFACES 

Figures 1(a) and (b) show our calculated band struc- 
tures of Ag(lll) and Au(lll) surfaces, respectively. 
Compared with corresponding experimental valueSfi^ the 
maximum binding energies of our surface bands are 
60 meV larger for Ag(lll) [Fig. 1(a)] and 100 meV 
smaller for Au(lll) [Fig. 1(b)], consistent with previous 
all-electron calculations The accuracy of the binding 
energy can be improved by careful description of de- 
cay of the wave functions into the vacuum,— but we 
already have good accuracy for the Rashba-type split- 
ting. Our surface bands show fco = 0.013 A^^ for 
the Rashba splitting in k space, defined by the offset 
by which the extremum of the surface-state dispersion 
shifts from the crossing point F and = 2.8 meV 
for the Rashba energy, defined by the energy difference 
between the crossing point and the maximum binding 
energy. When we fit our surface bands to an expres- 
sion, E± = ± afc -|- E'o, we obtain the electron ef- 
fective mass m* = 0.26 me, where mo is the electron 
bare mass, and the Rashba parameter a = 0.460 eV-A. 
These results are in good agreement with the reported 
experiments^ '^^i^^ and theoriesJ^ii^i^ In addition, if we 
shift the surface bands to match the maximum bind- 
ing energy with the experimental value, 487 mcV for 
Au(lll), the calculated Fermi wave vectors are 0.165 
and 0.194 A~^ for the inner and outer Fermi surfaces 
(FSs), respectively, which are also in good agreement 
with reported experimental and theoretical resultsi ^^i^^ 
These results show that our Kohn-Sham (KS) Hamilto- 
nian, HKsik) = Ho{k) + Hso{k), which is the sum of 
the scalar relativistic part, ifo(fc), and the SOI term, 
Hso{k), describes the Rashba-type splitting in Au(lll) 




FIG. 1. (Color online) Electronic band structures of (a) 
Ag(lll) and (b) Au(lll) surfaces around the 2D BZ center, 
f , calculated with the full SOI, Hso{k), and those of Au(lll) 
surface calculated (c) with only p-channel SOI, -ffgo^(fc), and 
(d) with only d-channel SOI, Hi^Z^ik). The Fermi energy is 
set to zero. The bulk-projected electronic states are shaded 
in cyan. 



surface accurately. 



IV. ORBITAL CONTRIBUTIONS TO THE 
RASHBA SPLITTING IN THE SURFACE 
STATES ON Au(lll) SURFACE 

In semiconductor devices, the relativistic response of 
electron spin to the electric field. Eh, yields the fc-linear 
Rashba Hamiltonian^ = aR(n xk)-(T, with the Pauli 
matrices a and api = efi^ E / (Amlc^) . These terms are 
derived from the microscopic SOI Hamiltonian 



Hi 



so 



:{VV xp)-a 



(1) 



and the nearly-free electron nature for 2D electron gases 
confined to surfaces or interfaces. While this model can 
explain the fc-linear Rashba splitting, it cannot explain 
the correct size of splitting^iii^ due to ignoring the effect 
of ion cores 1^ 
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To find the origin of tlie Rasfiba-type splitting, we con- 
sider orbital-angular-momentum (l) decomposition of the 
SOI term, i.e., Hsoik) ^ H'g^^{k) + H'g^^{k), and calcu- 
late the band structures with either i7gQ^(fc) or Hg^{k). 
This approach is possible because Hso{k) in our calcu- 
lation is expressed as a sum of ^-dependent Kleinman- 
Bylandcr-typc projectors. Figures 1(c) and (d) show that 
when only the p-channcl SOI, i.e., iJgQ^(fc), is switched 
on, the band splitting of the surface states is reduced 
greatly, but when only the d-channel SOI, i.e., H^{k), 
is switched on, the splitting is even slightly larger than 
that for the full Hso case. While m* in both cases are 
very close to the full Hso case, the Rashba parameter is 
ttp = 0.118 eV-A in the former case and ad = 0.564 eV-A 
in the latter case, giving the following relation 
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afuii 



ad 



(2) 



This shows that the Rashba- type splitting in Au(lll) 
surface originates from the d-channcl SOI while the p- 
channel SOI slightly reduces the splitting. 

In the above, we examined / decomposition of the 
SOI Hamiltonian. Now we consider orbital decompo- 
sition of the surface states on Au(lll) surface, i.e., 
V'nfc = 'ipnk,s+'4'nk,p + '4'nk,d, whcrc ipnk is the wavc func- 
tion of the nth surface band for the wave vector k and 
'>Pnk,s, tpnk,p, and 'ipnk,d arc s-, p-, and d-orbital parts 
of the state, respectively. With this decomposition, we 
calculate the expectation values {tpnk,i\Hsoik)\tpnk,j) of 
the full Hso for i, j = s, p, d to analyze orbital con- 
tributions to the Rashba- type splitting. We find that 
the d-d SOI energy, i.e., {ipnk,d\Hsoik)\'ip„k,d) , has oppo- 
site signs for the two bands, contributing most substan- 
tially to the band splitting. Then, the p-d SOI energy, 
i.e., {■4>nk,p\Hsoik)\tpnk,d) , contributes to the splitting by 
about 10 % of the d-d one. In contrast, the p-p and s- 
p SOI energies contribute weakly and negatively to the 
band splitting, i.e., they reduce the splitting slightly, and 
the s-s and s-d ones are almost zero. These results again 
show that the d-orbital character of the surface states 
determines the Rashba-typc splitting in Au(lll) surface 
although the surface states arc mainly of p-orbital one. 

To verify this picture for the difference in the Rashba- 
type splitting in Au(lll) and Ag(lll) surfaces, we cal- 
culate the band structures without SOI and estimate d- 
orbital characters in their surface statesj^i similarly to 
the previous worksJii^ As shown in Figs. 2(a) and (b), 
the Ag(lll) surface states have negligible d-orbital part, 
but the Au(lll) surface states have relatively significant 
d-orbital oneJ^ii^ This difference, due to the smaller en- 
ergy separation between Au bd and 6s levels than that 
between Ag Ad and 5s levelsj^ confirms the importance 
of d orbitals in the Rashba-type splitting Jiii^ 

As a double check, we consider shift of the bulk d-band 
energies, that is, we modify the KS Hamiltonian matrix 
by adding 5Ed to its diagonal elements for d orbitals after 
self-consistency, and then calculate the band structure. 
As shown in Table I, with no shift (6 Ed = 0), the Au(lll) 
surface bands arc split by A^p = 135 meV at the inner 
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FIG. 2. (Color online) Orbital characters of Au(lll) and 
Ag(lll) surfaces, (a, b) Band structures without SOI for (a) 
Ag(lll) and (b) Au(lll) surfaces. Both the color and the 
size of filled dots represent the norm of all d-orbital charac- 
ters (IV'nfeaiidl) '-'f sach wave function. The norm is defined 
using the pseudoatomic orbitals.— The bulk-projected elec- 
tronic states are shaded in cyan, (c, d) Orbital characters (c) 
without SOI and (d) with SOI for Au(lll) surface states along 
the positive kx axis. In (c), norms of s-, p^-, and p^-orbital 
characters (|^/'„e,J, and \'>p„k,p^) are in the left panel, 

and norms of d^z- and ci22-orbital characters (|i/'^j ^ \ and 
\^nk d 2 1) ^"^^ right panel. In (d), norms of d^z- and 

ci22-orbital characters (!?/'„£ ^ | and l-i/i^g d 2 D ''^ higher- 
and lower-energy surface states are in the left and right pan- 
els, respectively. The norm, IV'^g^li of an orbital character, 
is defined using the pseudoatomic orbitalsi^ 



Fermi wave vector, and the Rashba splitting in k space 
is fco = 0.013 A~^, as mentioned above, and the Ag(lll) 
surface bands have almost zero energy splitting (3 meV). 
For Au(lll), when SEd = —0.5 eV, the energy splitting 
Afcp decreases to 92 meV and fco decreases to 0.009 A~^ 
and when 5Ed =1.0 cV, A^p increases to 191 meV and ko 
increases to 0.019 A"^ For Ag(lll), when 5Ed = 5.0 eV, 
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Afep increases to 48 meV and fco increases to 0.006 A~^. 
These results show that the Rashba-type sphtting in no- 
ble metals indeed depends on d-band energies. We also 
note that 6Ed = 5.0 eV makes the shifted d-band energy 
in Ag close to the unshifted d-band energy in Au, but the 
splitting in Ag(lll) is still smaller than that in Au(lll) 
due to the difference in the intrinsic atomic SOI strength. 

Focusing back to Au(lll), Fig. 2(b) shows that the 
d-orbital part of the surface bands grows monotonously 
as k increases from F. This suggests that even the linear 
k dependence of the Rashba-typc splitting may be de- 
termined by /c-dependent size of the d-orbital character 
that exists without SOL Thus, with more analysis, we 
find that only s, pz, dxz, and (1^2 orbitals are significant 
for surface states along the axis, when the z axis is the 
positive direction normal to the top surface. As shown in 
Fig. 2(c), c?a;z-orbital character increases from zero, as k 
increases along the k^ axis, while s-, p^-, and c?22-orbital 
characters decrease slowly. Considering the symmetry, 
surface states along different k direction, i.e., along the ky 
axis, have dyz orbital significantly instead of dxz- When 
we use 



/,(0) 



for dxz- and (i;j2 -orbital parts 



and i^ld^, 

of a surface state ■(/'i^'' (not including spin) without SOI, 
respectively, we obtain that the SOI-induced energy split- 
ting Ak of the two surface bands is given by 



so('=)l^fc 



\k\. 



(0) , 



(3) 



where H^Q{k) is SOI between up- and down-spin states. 



Thus, the linearly fc-dependent ip^^l^ makes the Rashba- 
type splitting Afe proportional to As SOI splits the 
two surface bands, the d^^-orbital character increases in 
the higher-energy surface state while it decreases in the 
lower-energy one [Fig. 2(d)]. SOI redistributes d-orbital 
characters in the surface states but it does not mix the 
surface states with bulk states significantly. We note that 



TABLE I. Dependence of the Rashba splitting on the bulk 
d-band energy. Values of Akp and ko are calculated with 
the shift of the bulk d-band energy by SEd- For a better 
comparison, Ak;p refers to the splitting at the same k which 
is the inner Fermi wave vector in the case of SEd = 0. 
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even without SOI, the dxz and dz^ orbitals in the surface 
states along the kx axis form a nonzero y component of 
the orbital angular momentum, which may be related to 
a recent study of the Rashba splittingi^ 



V. SPIN AND ORBITAL ANGULAR 
MOMENTUM OF Au(lll) SURFACE STATES 

Another important feature of the Rashba-type split- 
ting is the helical spin structure. To find the role of d 
orbitals in the spin structure, we calculate expectation 
values of the spin angular momentum (SAM) operator. 
In addition, motivated by a recent work on the topo- 
logical insulator)^ we also calculate the orbital angular 
momentum (OAM) by projecting the surface-state wave 
functions onto a set of localized functions defined in a 
sphere centered at the top-surface Au atom and then 
evaluating expectation values of OAM operators around 
the atom using the projected wave functions.— A half 
of the nearest neighbor distance, 1.44 A, is used for the 
sphere radius, and real spherical harmonics multiplied by 
a constant radial part are used for the localized functions. 

As shown in Figs. 3(a) and (b), obtained SAM and 
OAM of the surface states are in-plane and helical around 
circular FSs. In our results, SAM is mainly from pz or- 
bital while OAM is mainly from d orbitals. When viewed 
from the vacuum side toward the top surface, SAM is 
clockwise at the inner FS and counter-clockwise at the 
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FIG. 3. (Color online) SAM and OAM of Au(lll) surface 
states at the Fermi energy, (a) SAM and (b) OAM with the 
full SOI. (c) SAM with only p-channel SOI, H's^\ (d) SAM 
with only d-channel SOI, Hg^. Black arrows represent SAM, 
and blue and red arrows do OAMs from p and d orbitals, 
respectively. For clarity, OAM from p orbital is doubled. 
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FIG. 4. (Color online) Effects of doping on the Rashba-type splitting on Au(lll) surface, (a) Squared amplitudes of surface- 
state wave functions at the same fcp, the Fermi wave vector of the inner FS in the undoped case, in 0.5 electron-doped, undoped, 
and 0.5 hole-doped cases, shown from left to right panels. The color scale is in units of per A^. The surface Au atom is marked 
with a dot at the center of each panel, (b) Band splitting A^p vs doping. All splittings are taken at the same fcp as in (a). 



outer FS [Fig. 3(a)], which is consistent with the reported 
experiments^ and calculations When viewed from the 
same side, 0AM from d orbitals is clockwise at both FSs 
while 0AM from pz orbital is counter-clockwise at the 
inner FS and clockwise at the outer FS [Fig. 3(b)]. To- 
tal 0AM is clockwise at both FSs and has smaller mag- 
nitude than SAM. These features for 0AM are also in 
good agreement with the recent experimental and the- 
oretical reports;^ but are different from the p-orbital- 
based Bi2Se3 topological insulator where 0AM has op- 
posite directions for the two surface bands and has larger 
amplitude than SAM.^ 

We also study the roles of p- and rf-channel SOIs in 
determining the spin direction. As shown in Figs. 3(c) 
and (d), our result shows that when only the p-channcl 
SOI, i.e., Hg^^{k), is switched on, the spin direction is 
opposite to the full Hgo case given in Fig. 3(a), but when 
only the d-channel SOI, i.e., Hg^{k), is switched on, 
the spin direction is the same with the full ffso case. 
This result shows again that the d-channel SOI is indeed 
decisive in the Rashba-type splitting while the p-channel 
SOI has minor and opposite effects. 



VI. ELECTRONIC PERTURBATIONS TO 
Au(lll) SURFACE 

As a way to change the Rashba-type splitting in the 
noble-metal surface, we consider electron and hole dop- 
ing, similarly to Xe and Na adsorption experimentsi^ We 
performed self-consistent calculations with different num- 
bers of electrons. In the undoped case, the surface-state 
wave function is asymmetric at the surface Au atom, 
with greater amplitude in the vacuum side than in the 
bulk side [Fig. 4(a)]. Compared with this, electron (hole) 
doping makes the wave function less (more) asymmetric 
[Fig. 4(a)], resulting from decrease (increase) of the dz2- 
orbital character in the wave function due to change in 
the degree of the inversion symmetry breaking, affected 
by change of the potential near the surface. As shown in 



Eq. ([3]), the Rashba splitting in noble- metal surfaces 
depends on the size of dz^ orbital in the surface-state 
wave functions. Thus, change in the o?-orbital characters 
changes the energy splitting. For Au (111) surface, we 
obtain 40 meV decrease for 0.5 electron doping per sur- 
face Au atom, and 10 meV increase for 0.5 hole doping 
[Fig. 4(b)]. These changes are ~100 times the Rashba 
effect in semiconductor devices at electric fields corre- 
sponding to the doping concentrations. 



VII. CONCLUSION 

In conclusion, we investigated the Rashba-type split- 
ting in the Shockley surface states on Au(lll) and 
Ag(lll) surfaces, based on first-principles calculations 
including SOI. We showed that although the surface 
states have predominantly p-orbital character and only 
small c?-orbital one, interaction between d-orbital parts 
via the d-channel SOI produces the Rashba-type split- 
ting. The small d^z ((ij,2)-orbital part of the surface 
states, present even without SOI, depends linearly on k so 
that the atomic SOI results in the linearly fc-dependent 
Rashba-type splitting. The d-channel SOI determines 
the helical spin structure of the surface states although 
the spins are mainly from p orbitals. The p-channel SOI 
has minor and opposite effects in Au(lll) surface. Thus, 
the role of d orbitals in Au(lll) surface has some sim- 
ilarity to graphene where small 3d orbitals contribute 
dominantly to the SOTinduced energy gap^ but shows 
difference from p-orbital-based Bi2Se3-type topological 
insulators)^ As a way of varying the splitting size in 
Au(lll) surface, we discussed effects of electron and hole 
doping. 
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